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Abstract 

In this paper, influence of crowding by inert particles on the geminate reaction kinetics is theo- 
retically investigated. Time evolution equations for the survival probability of a geminate pair are 
derived from the master equation taking into account the correlation among all diffusing particles 
* , and the results are compared with those obtained by Monte-Carlo simulations. In general, excluded 
£^ \ volume interactions by the inert particles slow down the diffusive motion of reactants. However, when 
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the initial concentration of the inert particles is uniform and high, we show that additional influence 
of interference between reaction and correlated diffusion accelerates the transient decay of the survival 
probability in the diffusion controlled limit. We also study the escape probability for a non- uniform 
initial distribution of the inert particles by taking the continuous limit in space. We show that reaction 
yield is increased when the reaction proceeds in the presence of a positive density gradient of the inert 
particles which inhibits the escape of reactants. The effect can be interpreted as a cage effect. 
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I. INTRODUCTION 



The importance of molecular crowding on chemical reactions has attracted great attention 
in connection with biochemical reactions in living cells. jl-7| Living cells contain a high volume 
fraction of macromolecules, in addition to reactants. Although these macromolecules are not 
reactive, the excluded volume interactions between reactants and macromolecules significantly 
affect transport properties of reactants, and therefore biochemical reactions. 

In this paper, we consider a fundamental reaction process called geminate reaction, which is 
observed in many systems including those encountered in biology. In geminate reaction, a pair of 
reactants is generated simultaneously and subsequently diffuse and react when they encounter. 
Geminate reactions are influenced by spatial diffusion of a pair and the intrinsic recombination 
rates. The influence of many body interactions between inert species and reactants on the 
geminate reaction kinetics can be very complicated and difficult to treat theoretically. The 
simplest model could be to assume that the reactants and the inert particles have the same 
size. Even under such simplification, the many body nature of the problem remains since the 
migration of reactive species correlates with the time dependent positions of inert species; the 
problem is still difficult to solve analytically. 

In order to retain many body nature in the simplest situation, we study geminate reaction 
between a static species and a diffusive species on a lattice. Reaction takes place according to 
the distance between one of the pair of reactants at the origin and the other. Inert particles 
perform random walks on a lattice. The transition to neighboring lattice sites is constrained by 
prohibiting the double occupancy; each lattice site can be occupied at most by a single diffusive 
particle regardless of whether it is reactive or inert. Particles are assumed to move randomly 
on vacancy sites of a lattice. However, the problem is still hard to solve analytically without 
approximations. Therefore, we perform Monte-Carlo simulations to evaluate approximations 
and to elucidate effects which cannot be studied analytically. In order to facilitate comparisons 
between theoretical results and those of simulations, the problem is further simplified; the origin 
is also allowed to be occupied at most by a single diffusive particle regardless of whether it is 
reactive or inert and reaction takes place according to the intrinsic reaction rate when the origin 



2 



is occupied by a reactant. 

The excluded volume interactions were theoretically treated by Nakazato and Kitahara in 
tracer diffusion on a lattice. js| Nakazato-Kitahara's formula of tracer diffusion constant inter- 
polates between low and high concentrations of host particles and its accuracy is confirmed by 
comparison to the results of numerical simulations. |8l-lll| 

For target reactions where a static reactive particle (target) is surrounded by many reactive 
counterparts (quenchers), Nakazato-Kitahara's theory was successfully applied to calculate the 
survival probability of a target with a constraint of prohibited double occupancy of diffusing 



reactants. 
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13] It turned out that the decay of the target survival probability is accelerated 



by prohibiting the double occupancy. 
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other numerical and analytical approaches. 



Simi 



ar acceleration of the decay was obtained by 
14H25| The acceleration of the decay is understood 



by noticing that the number of sites occupied by mobile reactants is generally l arge r at any time 
under the constraint of prohibited double occupancy at each lattice site. [12, ll3[ Accordingly, 
the probability of reaction between the target and a quencher is higher at any time when 
multiple occupancy is not allowed. 

Contrary to the target reaction, only a pair of reactants should be considered for geminate 
reaction. In other words, the number of sites occupied by reactants is not affected by prohibiting 
the double occupancy. However, the site blocking effects among diffusing particles influence 
the kinetics of geminate reaction through different mechanism from that in the case of target 
reactions. First, crowding of inert particles slows down diffusion of reactive particles and retards 
the reaction between the pair. Indeed, numerical simulations show that the reaction between 
a pair proceeds slowly by the crowding of inert particles. 26] 

In this paper, we study more comprehensively the reaction of a pair of reactants performing 
diffusion under the constraint of prohibited double occupancy in the presence of inert particles. 
On the basis of results derived, the mechanism of site blocking effects by inert particles in 
the geminate reaction is investigated in detail when the initial distribution of inert particles 
is uniform. We show that reactions are influenced through excluded volume interactions not 
only through slowing down of diffusion of reactants but an interference between reaction and 
correlated diffusion in the diffusion-controlled limit. We also study the influence of inhomoge- 
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neous initial distributions of inert particles by taking the continuum limit in space. We show 
that the overall reaction yield is increased (decreased) from that assuming homogeneous initial 
distribution of inert particles by a positive (negative) density gradient of the inert particles. 

In Sec. II we formulate the problem in the case when the initial distribution of inert particles 
is uniform and derive the solution within the mean field approximation. In Sec. Ill higher 
order corrections to the mean field results are presented. The details of calculation are given 
in Appendix El In Sec. IV, we compare the analytical results to simulation results. In Sec. V, 
the influence of non-uniform initial distribution of inert particles is investigated. The escape 
probability is derived for this case by using the continuous limit derived in Appendix [B] Section 
VI is devoted for conclusions. 



II. GEMINATE PAIR REACTION UNDER THE PRESENCE OF INERT PARTI- 
CLES 



For simplicity, we formulate the problem on a lattice where a reactive particle and inert 
particles perform random walks. The particles are assumed to move randomly on the vacant 
sites of a lattice. One of the reactants of the pair does not move and its position is taken as 
the origin of the coordinate system. The reactive particle undergoes reaction according to the 
distance from the origin, r. We denote the intrinsic reaction rate by k(r). 

The tracer-diffusion in concentrated lattices was studied by Nakazato and Kitahara in the 

n 

absence of reaction. |8| The diffusion of the tagged particle in the presence of site blocking 
by other particles has been studied. jsil] Following them, we introduce ket vectors to show 
occupancy of a site by diffusing particles. The ket vector \ f, •) denotes the occupation of site r 
by a reactive particle, the ket vector \r, o) denotes the occupation of site r by a inert particle, 
and \r,<f>) represents that site f is empty. The conditional probability of finding inert particles 
a t /2j " " " j Jn) and the reactant at n at time t when the initial configuration of inert particles 
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is (ii, ii-, • ■ ■ , in) and that of the reactant is fh is written as, 



iV 



M 



P(fi,f2r--,fN,n,t;i 1 j 2 ,---j N ,m) = lY[{f£, \)[ {hA\ (n,«|exp(£t) 



N 



M 



K-) ]Il4o) 



=JV+1 



where iV and M denote the numbers of inert particles and lattice sites, respectively. L is given 
by the sum of the term describing diffusion L w and that describing reaction L rc , C = L w + L rc . 
L w is explicitly expressed as, |s-10| 



|r n ,«)(r n ,»| • \ f m , (f)) (f m , (f>\) + 



(n.m) 



(n,m) 



(2) 



where the sum is taken over all nearest neighbor pairs of the accessible lattice sites by the 
diffusing particles. 7b is given by 7# = Tb/ (2d), where is the jump frequency of a reactive 
particle and d denotes the lattice dimension. Similarly, we define ^ w = T w /(2d) where T w is 
the jump frequency of inert particles. L rc describes the reaction from an occupied site f n with 
the rate k(r n ), |27H29| 



L rc — ~ } k (r n ) \r n , *){r n , 



(3) 



n=l 



The conditional probability, Pn {n, t\fh, 0), that the reactant is at site n at time t when it 
was initially at m under the assumption of random initial occupation of inert particles is 
obtained from Eq. (JTj) by multiplying 1/ (mCn) and summing over all possible initial and final 
configurations of the inert particles. By defining the characteristic function by, 

M 

,t\m,0;x) = ^2P N (n,t\m,0)x N , (4) 



in, 



it can be expressed as, 



[n, t\m, 0; x) 



N=0 



[l+x 



M 



mCn 



-g [n, t\m, 0; x) 



(5) 
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where we define, 

g (ft, t\m, 0; x) = ({0} | (ft, • exp (L(9)t) | {</>}) \m, •), (6) 



L(6>) = exp (— ^S 1 ) £exp (^S 1 ), S 1 = ^^ =1 0| — °|), and x = tan 2 6>. It is 

convenient to introduce abbreviations, 



' M 



where M' denotes that the site fj is excluded in the product. 

The inverse transformation is given by applying the Cauchy's integral theorem, 

I f 1 

Pat (n,t\m, 0) = — / dx-^^(f)(n,t\m,0;x), (9) 

where the path of integration encircles the origin on the complex plane. 

In the thermodynamic limit of M — > oo with c = N/M being constant, the right hand side 



of Eq. (jH]) can be calculated by applying the saddle point method, |8HlO| 

P (n, t\m, 0) = g (ft, t\m, 0; c/(l - c)) . (10) 
The survival probability of a pair at time t whose initial separation is given by m is defined by, 

S{m,t) = ^2P{n,t\m,0). (11) 

From Eqs. and (TIT)]) the Laplace transform of the survival probability, S (m, s) = 

J °° dtexp(—st)S (m,i), is expressed as, 

S (m, s) = \(n, -I ly— A M)K •>, (12) 

~^ s-L{6*) 



where tan#* = y c/(l — c). L{6*) can be expressed by the sum of the term describing diffusion 
and that describing reaction, L{9*) = L rw + L TC . Even after the transformation the term 
describing reaction is not changed, L TC = L rc , while L rw is given by the sum, L rw = L Tw0 + L TW i, 

6 



w 

[8 



lere L rw o describes the transition under the conservation constraint of the number of particles, 
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L rw0 =(1 -c)7b ^ (\r n , •) (r n , <fi\ • \r m , (f>) (r m , •] - K,»)(r n ,«| • \r m ,<f>){r m ,<f>\) + 

(n,m) 

ciB ^2 (\^n, m ){r n ,°\ ' \r m , °) (r m , »\ - \ f n , •)(f n , •! • \r m ,o){f m ,o\) + 

(n,m) 

lw (|f„,o)(r n ,0| ■ \ f m ,(j)){f m ,o\ - \ f n ,o)(f n ,o\ ■ \r m ,(j)){f m ,(p\) , (13) 

(n,m) 

and L rw i describes the transition where the number of particles is not conserved, js-10] 



L rw i =a/c(1 - c) 7^ ^ (|f n , •)(f„, •] • \ f m ,(f)){f m , o| + 

(n,m.) 

r n ,«){r„,«| • |r m ,o)(f m ,0| - |f n , •)(r„,0| ■ |r m ,o)(r r , 
|r n , •)(f„,o| • |f m ,0)(r m , «|) . 



(14) 



By introducing the identity, 



s-L( 



- ( 1 + 1 L(9* 

s V s-L(9*) 



(15) 



Eq, (TT2j) can be rewritten as, 

sS (m, S ) = 1 + Ift H ^7-7^(^)1 M>K •)• 



(16) 



In the lowest order approximation, the perturbation term, L rw i, is ignored in the numerator 
of Eq. f)16p and we obtain, 

1 



sS(m,s) = l + J2m 



s-L(6*] 



L rw0 + L rc | {0})|m, •). 



(17) 



By using Eq. ( fl2i) and the fact that the number of particles is conserved for both L rw o and L rc . 
Eq. (HZD leads to 

sS (m, s) — 1 = 7_b(1 — c) ^S 1 + bj, — (m, s) j — k (m) S* (m, s) , (18) 



where fh + bj denotes a nearest neighbor of the site fh and the sum is taken over all nearest 
neighbor sites. By the inverse Laplace transformation, the equation for the survival probability 
at time t of a pair with initial separation fh is obtained, 

—S (fh, t) = 7B (1 -c)J2(s(m+ bj, tj - S (rh, *)J - k (rh) S (fh, t) . (19) 

3 

In the lowest order approximation, the site blocking effects by inert particles reduces the tran- 
sition rate. The transition rate is reduced since jump to a neighboring site is allowed only 
when the neighboring site is empty. The vacant probability is 1 — c in the mean field picture. 
Equation ( I19j) is a mean-field result in the sense that the reduction factor is given by 1 — c. 
The transition rate of the reactant particle decreases linearly with increasing the concentration 
of inert particles. 

For localized reactions, k (fh) = k^S^ the general solution after the Laplace transformation 
is obtained as, 

or-' \ 1 ( i G (fh,s)k 

S(m,s) = - 1-- A t7 . - , (20) 



s 



l + G (0, S )fc 



where the Green's function, 



G (j,s)= 1 Ms) U(j,s), (21) 
s 



is given in terms of the lattice Green's function, 

1 r /-TV exp ( — ik ■ j 
U(j, s) = ■■■ d d k > (22) 

u ' ; (^yi ]_ v i-^ B ( s )\(k) y } 

where ^s(s) is given by iPb(s) = r B (l — c)/(s + T B (1 — c)), the structure factor is defined by 
X(k) = Y^j=i cos (k ' bj/bj and b denotes the lattice spacing. 

The recombination probability of a particle starting from fh, k (rh) = 1 — lim^oo S (rh, t), is 
obtained as 

U(fh,0) , , 

«(^)= r fi c) • (23) 

kn 
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Note that U(m, 0) for any fh is independent of the concentration of the inert particles, c, in 
the mean-field result. In the limit of perfectly absorbing boundary condition (ko —> oo), the 
recombination probability is independent of the concentration of inert particles. For partially 
absorbing boundary conditions the recombination probability increases by increasing the con- 
centration of the inert particles. The escape probability, liim-^ S (fh, t), which is defined as 
the probability of a pair with initial separation fh surviving at infinite time, is given by, 

ip (m) = 1 - k (fh) . (24) 

We have derived the simplest results on the survival probability of a geminate pair by ignoring 
correlations higher than the two-point correlation between the initial position and the position 
at an arbitrary time. In the reaction-diffusion equation thus derived, the presence of the inert 
particles only reduces the diffusion coefficient linearly with increasing the concentration of the 
inert particles and the diffusion and the reaction do not interfere. In the subsequent section, 
we show that the diffusion and the reaction interfere in the presence of inert particles if we 
consider higher order correlations. 



III. CORRECTION TO THE MEAN FIELD EQUATION 



If we ignore correlations higher than two-point correlations, the Bardeen-Herring back cor- 
relation is not taken into account. 3l| The Bardeen-Herring back correlation takes place when 
the reactant hops to a vacant site, leaving the previous occupied site vacant; after the hopping 
the transition probability of the reactant back to the previously occupied site is higher than 
other sites. The velocity autocorrelation function in a lattice gas shows a long time tail with 
a negative value due to the Bardeen-Herring back correlation, Ill] Suppose that a reactant 
occupies a reactive site after a hopping. The rate of hopping back to the previous site competes 
with that of reaction. In this way, the reaction interferes the correlated diffusion. Interference 
means that the reaction process and the diffusion are not statistically independent. In this 
section, we study the interference between the reaction and the correlated diffusion by taking 
into account the Bardeen-Herring back correlation. As in the previous section, we assume the 
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initial uniform distribution for the inert particles. 

The effect of the interference between the reaction and the correlated diffusions can be 
calculated as the correction to the simple diffusion-reaction equation, Eq. (!T7|) . The exact 
relation, Eq. ffTB"]) . can be rewritten as, 

sS (m, s) = 1 + | (ft, «| 1— (L rw0 + L rc ) | {0}> \m, •) + R(m, s), (25) 

^ s - L{6*) V / 

where R(fh, s) represents the correction to Eq. (ITT)) and is given by, 

R(m, s) = J2({<P}\ (n, • I ^T^wi | W) |m, .) . (26) 

By noticing L rwl = L(0*) — L rw0 — L rc , we can prove the operator identity, 

I— - = l—L^ ^ =r- + ^ ^. (27) 

s — L[9*) s — L[8*) s — L tw q — L rc s — L tw q — L TC 

L rw o and L rc conserve the number of • in the bra and ket notations, while L rw i does not, so we 
have 

({0} |(n, .| ^ ^L rwl | {0}}|m, •) = 0. (28) 

If we substitute Eq. ( 1271) into Eq. ( 1261) and use Eq. ( 1281) . Eq. ( f26|) can be expressed as, 

R(m,s) = ^(W|(n>| FT^^rwi r ^ ^A-wil (29) 

H S — L{U ) S — L Tw0 — L rc 

= Y]S (n,s) ({4>} \ (n,»\L rw i ^- ~^L Twl \ {4>})\rh, •), (30) 

s — L 



rwO -^rc 



where the definition of S (n, s) given by Eq. ( 1T21) is substituted. By introducing the explicit 
expression of L rwi given by Eq. (EH]), we obtain, 

({0} 'l-^rwl pi ^^rwl| {0})|m, •) = 

s — L rw o — L IC 



lie 



(1 — c) [g n + b r o \rh», fa + 6 9 o, s j — G + fe r «, n o |m«, m + fe g o, s 

r q 

—G ^n«, n + 6 r o |m + mo, s J + G (n + 6 r «, n o \ fh + mo, s^j . (31) 
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Here we define the four-point correlation function as, 

G(fi;rj o \r k ;r e o,s) = ({(f)} \(fi, •|(r i , o| =— =— | {0})|r fc , m)\f e , o) (32) 

s — L rw o — L rc 

using the abbreviation, 

(m" \ 
n(rl,0|j (rl,.|(rS,o|, (33) 

where M" denotes that the site r*j and the site r} are excluded in the product. Eqs. (!33l) - (l34l) 
represent the state that all sites are vacant except the site r*j occupied by the reactant and the 
site fj occupied by an inert particle. By substituting Eqs. (13"Uj) and (l3"Tj) in Eq. (1231) . we have, 

sS (m, s) — 1 = 7s(l — c) -^(^j &d s ) f'S' (™ + &rj s j — S 1 (ft, s)J — (m) (m, s) , 

n r 

(35) 

where the kernel F c (n,m,b r , s) is given by 
F c (n, m, 6 r , s) = 5^ - 7 B c ^ n + b r o \ fn», m + 6 g o, sj 



— G (ft*, n + 6 r o |m + b q », mo, . (36) 



By the inverse Laplace transformation of Eq. ( 1351) . the survival probability is shown to satisfy 
the diffusion-reaction equation in which the diffusion term is expressed by the time convolution 
with the nonlocal kernel, F c (n,m,b r ,t). 

The time convolution represents the memory effect originating from the correlation between 
the mobile reactant and the inert particles. Reactant motion is correlated with the time de- 
pendent arrangements of the inert particles through prohibited double occupancy of the lattice 
sites. In particular, the site occupied by the reactant becomes empty just after the hopping 
of the reactant and the chance of back transition to the previously occupied site is high. The 
back transition probability of the reactant decreases as time proceeds because the empty site 
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generated by the hop of the reactant to a vacant site may be occupied by another inert particle. 
The time dependence of back-jump correlation is the origin of the memory effect. 

In principle, the back-jump correlation competes with reaction. Suppose that the reactant 
hops to the reactive site. The probability of jump back to the previously occupied site decreases 
as the reaction rate increases. Since Eq. (|32|) includes the operator describing reaction, L rc , 
the diffusion memory kernel, F c (n,rn,b r , s), depends on the reaction rate. In order to obtain 
F c (n, m, 6 r , s) we need to solve an equation for G (n», n + b r o |m«, ro, s j . 

In Eq. (1A2|) of Appendix |AJ we show that the equation for G (^n»,n + b r o |m«,ro, s^j in- 
cludes the reactive sink term. The interference between reaction and correlated diffusion is 
taken into account by the four-point correlation function. In the simplest theory given by the 
two-point function, Eq. (fl9]) . the interference between reaction and correlated diffusion is not 
taken into account. 

When the reactive sink strength changes according to the distance from the origin, the dif- 
fusion term given in terms of the four-point correlation function depends on the distance from 
the origin accordingly. In addition, the presence of the inert particles gives rise to correlation 
over distances as a result of the excluded volume interactions between the inert particles and 
the reactant. Interference between reaction and the correlated diffusion beaks down the trans- 
lational invariance as shown in Eq. f ]A2j) of Appendix |A] and the resultant equation is hard to 
solve. In the next section, we use numerical simulations to study the interference effect. 

When we ignore the interference between reaction and correlated diffusion, the translation 
invariance is satisfied for the four-point correlation functions. Under the translational invari- 
ance, G (n», n + b r o \fn», ro, s^j depends only on relative vectors and satisfies, 

G (n», n + b r o \rh», ro, = G (n'», n' + b r o |0«, r'o, sj , (37) 

where n' = n — fh and r' = r — m. Since the number of independent variables is reduced, it is 
convenient to introduce a new notation, 

G [T) (n'm,b r o\f'o,s) = G (n';n' + b r o |0«,r'o,s) . (38) 

In this case, the time evolution equation for the survival probability is expressed after the 
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spatial Fourier transform as, 

(k, tj = J dtiM(k, t-h)S (k, h) - ^ exp (ik ■ fnj k (rh) S (rh, t) . (39) 

In the Laplace domain, M(k, s) can be regarded as a self-energy or memory function and is 
expressed as, 

M(k, s) = -r fl (l - c) (l - A(fc)) + 7 Ic(l - c)f c (k, s) , (40) 

where the correlations among the inert particles and the reactant as a result of the excluded 
volume interactions are included in, 

f c (k, s) = i 1 - ex P • ^)) ^ (X &r |6„ s) (l - exp (ifc • 6,)) , (41) 

q r 

G (jfe, 6 r |f, s) = ^ exp (iifc • /) G (T) (7, 6 r |f, s) . (42) 

The same form of memory function expressed in terms of the four-point correlation function 
was derived by a different method. [32] The equation for G (k, b r \r, s\ is explicitly shown in 
Eq. (JSTJ) of Appendix El 



In the limit of small wavelength, k — > 0, Eq. (139]) can be expressed after the inverse Laplace 
transformation as, 

(rh, t) = D B (1 - c) J2 f dtj c (t - h) V 2 S (rh, h) -k(rh)S (rh, t) , (43) 

at . j 

where the diffusion constant is defined by, Db = b 2r jB, the correlation factor in the Laplace 
domain is given for the hypercubic lattices by, 

f c (s) = l- lBC J2 [G (a) (br \K, s) - G (a) (br | - b r , s)] , (44) 

r 

and the initial condition is S (rh, t — 0) = 1. G^ (b r . \b g , s\ is defined by, 

G (a) (b r K s) = K", + K, o\ * 1 {</>}> |m, .)|m + 6„ o), (45) 



rwO 
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which is independent of the choice of n as shown in Eq. (IA14I) of Appendix The equation for 
(b r \b q , s^j is known and has been studied to obtain the tracer-diffusion coefficient. 8-1(3] 
It is shown in Eq. (1A14I) of Appendix |A] Its solution is known and the correlation factor can 
be expressed as, 

« s) = : . ^A-acr (46 > 

1 — M s ) r-i r 

V ; 1 w + 1b(1-c) 



1 f r , d ? 2sin 2 fc 1 



where is given by, 

u(s) = ~, ~~t~7 / ■ ■ • / d d k z^L^l _ (47) 

and 7t = 7«, + 7b(1 — c). In the original derivation of the correlation factor f c (s) and the 
tracer-diffusion coefficient by Nakazato and Kitahara, a projection operator method is applied. 
Here, a reaction-diffusion equation for the survival probability is derived directly without using 
the projection operator method. 

In order to take into account the memory effect on the transient decay of the survival 
probability with reasonable simplicity, Eqs. (|20|) - (|22|) are used with 

Us) = • (48) 

s + r B f c ( s )(i-c) 

In this approximation, the lattice Green function valid for finite wave length is used together 
with the expression of / c (s) derived in the limit of k — > 0. In the subsequent section, we show 
by comparison with simulation results that the approximation gives reasonable results as long 
as the interference between reaction and the Bardeen-Herring back correlation is absent. 

The equation can be further simplified by ignoring the memory in the diffusion kernel. In 
this approximation, the survival probability and the recombination probability for localized re- 
actions can be calculated respectively from Eq. ( 120]) and Eq. ( 123]) by introducing the correlation 
factor into the hopping frequency, 

1b -»■ JbU (49) 
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where the correlation factor f c = f c (s = 0) is given by 



fc 



1 — /i 



(50) 



lw + JbQ- - 3c) ' 



and /i = /i(s = 0) is known for some lattices. When the hopping frequency are the same for 
the inert particles and the reactant, 7^ = 75, the value of \i is 0.20984 and 1 — (2/tt) = 0.363 
for the cubic and the square lattice, respectively. 

The summary of this section is as follows. We have studied the influence of back-jump 
correlations on the survival probability of a geminate pair when the initial distribution of inert 
particles is uniform. When the mobile reactant hops to a vacant site, the reactant tends to jump 
back to its previously occupied empty site (the Bardeen-Herring back correlation). In this way, 
the reactant motion is highly correlated with the time dependent arrangements of the inert 
particles. The back-jump correlations interfere with reaction. In principle, the interference 
can be taken into account by Eqs. (135]) and fl36|) with Eq. (IA2I) . However, in practice these 
equations are hard to solve. If the interference is ignored, the influence of back-jump correlations 
is taken into account by Eqs. fl39l) - fT4^1) with Eq. (IA12I) . By introducing further simplification 
of ignoring the memory effect, we obtain Eq. (TT9~|) with substitution given by Eq. (l49l) . In the 
next section, these results will be compared with those obtained by numerical simulations. 

IV. COMPARISON TO SIMULATION RESULTS 
A. Simulation method 

In order to see the interference of reaction with the correlated diffusion, we perform Monte- 
Carlo numerical simulations. We numerically obtain the probability of geminate reaction in 
the presence of site-blocking effects using a kinetic Monte Carlo method. The simulation is 
carried out on the simple cubic lattice. One reactant is placed at the lattice site (0, 0, 0) 
and assumed to be immobile. The other reactant is initially placed at (j, 0,0), where j is an 
integer and the lattice constant is unity. Inert particles are randomly generated at lattice sites 
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within the box (0, L — l) 3 , where L is the box length. The number of inert particles, N, is 
related to their concentration by c = N/L 3 . Each lattice site may accommodate only one 
inert particle or the mobile reactant. We assume that the inert particles belonging to the box 
(0, L — l) 3 are periodically replicated in three dimensions, so that the simulation volume is 
effectively unlimited. What should be noted is that the spatial periodicity is assumed only for 
the distribution of inert particles, and not for the reactants themselves. During the simulation, 
both the inert particles and the mobile reactant may perform hops to neighboring lattice sites. 
A hop is allowed only when the destination site is not occupied by another inert particle or the 
mobile reactant. However, it is allowed for both types of simulated particles to jump to (0, 0, 0). 
If the mobile reactant is staying at (0, 0, 0), its reaction with the other reactant is possible. The 
procedure of selecting the event that actually takes place at a given simulation step is as follows. 
First, we determine all possible hops for the mobile reactant and inert particles. Denote the 
numbers of such hops as K and K m , respectively. If the mobile reactant is staying at a site other 
than (0, 0, 0), the total rate of all possible events is calculated as k to t = K^B+K- in ^ w . Otherwise, 
the total rate includes the rate of reaction and is calculated as k tot = K'js + -KinT™ + ^o- Now, 
we determine which event will actually take place. This is decided at random, taking the 
ratio of the rate of each possible event to the total rate k to t as the event probability. The 
above procedure is repeated until either a reaction occurs or the mobile reactant separates to 
a large distance r max from (0,0,0). By repeating the simulation for a large number (at least 
2 x 10 4 ) of independent runs, we can obtain the reaction probability. The accuracy of the 
simulation results depends on two parameters: L and r max . They should be taken as large as 
possible within the practical limits imposed by the available computational time (the demand 
on computer time is especially high at large concentrations of inert particles). In the production 
runs of the simulation, we assumed L = 10 and r max = 30. From test calculations carried out 
also for other values of these parameters, we found no significant effect of L on the obtained 
results. However, a weak dependence of the reaction probability on the value of r max could be 
observed. For example, the reaction probability obtained for j = 3 and c = with r max = 60 
was about 2% higher than that calculated with r max = 30. 



1(3 



B. Simulation results 



We investigate quantitatively the effects of the factors ignored in deriving simple result, Eq. 
(j!9p . by comparison to the more rigorous theoretical results and the simulation results. 

One of the factors ignored is the Bardeen-Herring back correlation. The Bardeen-Herring 
back correlation is described by the four-point correlation function given in Appendix |X] The 
Bardeen-Herring back correlation is taken into account fully by Eqs. (j39j) - (l42j) and partly by 
Eq. flU]) with Eq.fllHD. Equation §M) is obtained from Eqs. (EH])-(S2D by taking the limit of 
small wave length, k — > and ignoring the memory effect. Equation (Tl9l) with Eq. P9|) is much 
simpler than Eqs. ( )39|) -( l42|) . The numerical way to solve Eqs. (I39|) - (j42|) with the additional set 
of equations are given below Eq. (IA12j) of Appendix lAl 

Another factor ignored is the effect of the interference between reaction and the Bardeen- 
Herring back correlation. The interference is taken into account in the results of numerical 
simulations but is ignored in any theoretical results including the most sophisticated one given 
by the solution of Eqs. §9§- @2}- 

The recombination probability in the completely diffusion controlled limit, ko — > oo, is shown 
in Fig. [TJ The simulation results are compared to the numerical solutions of Eqs. (T39]) - (142]) . 
In the theoretical results the influence of reaction on the four-point correlation function is 
ignored. As shown in Fig. (TJ the influence of reaction on the four-point correlation function is 
relatively small for all concentrations of inert particles as long as overall yield (recombination 
probability) is concerned. However, Fig. [I] shows small deviation at high concentration of inert 
particles for the initial separation of j = 1. The deviation is not an effect of statistical error of 
the simulation. We will discuss this point later when we study transient decay of the survival 
probability. 

In the completely diffusion-controlled limit, ko — > oo, the recombination probability obtained 
from Eq. ( )23l) with the substitution given by Eq. (H9l) is independent of the concentration of the 
inert particles. This result is not a rigorous relation. It is obtained by the oversimplification of 
taking the limit of small wave length, k — > 0, and ignoring interference between reaction and the 
Bardeen-Herring back correlation. However, the difference between the result of the simplified 



17 



0.0 1 ' ■ ■ ' 1 ' ' ■ ' 1 

0.0 0.5 1.0 

c 

FIG. 1: Recombination probability of a geminate pair against concentration of inert particles, c, for 
ko — > oo (In the simulation, fco/7s = 10 6 ) and r y w /^/B = 1- j indicates the initial separation of the 
geminate pair. Dots with error bars indicate the simulation results. Squares represent the numerical 
solutions of Eqs. (|5g ]) -(]35 ]) with Eq. (|A1~2|) . Dashed lines indicate the mean field results of Eq. ()23[) . 
The results of Eq. (|23[) with the substitution given by Eq. (|49p practically coincide with those of Eq. 

(ESD. 

equation, Eq. with Eq. P^l) . and that obtained without taking the limit of long- wave 

length, Eqs. ( |39]) - P2j) . is also very small. Figured] indicates that the simplified approach which 
leads to Eq. ( )23|) with Eq. (|49|) is justified for the calculation of the recombination probability 
in the diffusion-controlled limit. 

The recombination probability in the case of finite reactivity is shown in Fig. |2j The 
influence of reaction on the four-point correlation function can be seen as the difference between 
the simulation results and the numerical solutions of Eqs. ( I39p - (H2|) . In the reaction-controlled 
limit, the difference is very small for all concentrations of the inert particles regardless of the 
initial distance of a geminate pair. The difference between the simplified results obtained Eq. 
( I23|) with Eq. (1491) and the solutions of Eqs. (!39|) -( l42|) is again negligibly small. 

We also compare transient decay of the survival probability obtained from simulations with 
theoretical results of Eqs. fl20"|) - (l22|) with Eq. (H8"j) to study the effect of interference between 
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FIG. 2: Recombination probability of a geminate pair against the concentration of the inert particles, 
c, when r y w /'fB = 1- (a) ^o/tb = 1-0 and (b) fco/TB = 0.5. j indicates the initial separation of the 
geminate pair. Dots with error bars indicate the simulation results. Squares represent the numerical 
solutions of Eqs. (I39p -(|42p with Eq. (|A12j) . The results of Eq. (|23|) with the substitution given by 
Eq. (|49p are indicated by the solid lines. Dashed lines indicate the mean field results of Eq. (|23p . 



reaction and the Bardeen- Herring back correlation. Figure [3] (a) shows that the theoretical 
results reproduce simulation results in the reaction-controlled limit (fco/7s = 1-0) over the 
whole time range regardless of the concentration of inert particles. As we have theoretically 
shown in the previous section, the interference between reaction and the Bardeen-Herring back 
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FIG. 3: Survival probability as a function of 7#t for the initial separation of j = 1. (a) ko/^B = 1-0 
and (b) /co/tb = 10 6 . Concentrations of inert particles are shown in figures. Symbols represent the 
simulation results. Lines indicate the theoretical results of Eqs. (|20p - (|22p with Eq. (|48p . 

correlation changes the diffusion term through 4-point correlation functions, leaving the reaction 
term unaltered. The interference could be dominated in the diffusion-controlled limit but should 
be small in the reaction-controlled limit. Figure |3] (a) indicates that the approximation used 
for finite k in Eqs. (120|) - (|22|) and Eq. (148]) is valid as long as interference between reaction and 
the Bardeen- Herring back correlation can be ignored. Figure |3] (b) shows that in the diffusion- 
controlled limit the theoretical results agree with the simulation results in the absence of inert 
particles. However, when the concentration of inert particles is high (c > 0.5), the survival 
probability obtained from simulation decays faster compared to the theoretical results. The 
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acceleration of the decay should be attributed to the influence of interference between reaction 
and the Bardeen- Herring back correlation. 

As the conclusion of this section, we point out that the interference between reaction and 
the Bardeen- Herring back correlation is most pronounced in the transient decay of the survival 
probability in the diffusion-controlled limit at high concentration of inert particles. As long 
as the reaction yield is concerned, the mean field results given by Eq. (fT§l) reproduce the 
simulation results regardless of reaction strength and the concentration of the inert particles 
when the density of the inert particles is uniform. The result can be improved by introducing 
substitution given by Eq. fH9|) into Eq. (I19p to take into account the Bardeen-Herring back 
correlation. Generalization of Eq. (TT9]) to the case of continuous diffusion and non-uniform 
distribution of inert particles is shown in the subsequent sections. 



V. INHOMOGENEOUS DISTRIBUTION OF INERT PARTICLES 



So far, we have assumed the homogeneous distribution of inert particles. Recently, the influ- 
ence of inhomogeneous distributions of inert particles is taken into account to study catalytic 
surface reactions, in particular focusing on reaction front structures. |33l-l38| In this section, 
we study geminate reactions under the inhomogeneous distribution of inert particles. Since we 
are not able to solve the lattice model for nonuniform initial distribution of inert particles with 
the same rigor as that under uniform distribution of inert particles, we study the results in 
the continuous limit by using the mean field approximation. We consider the pair distribution 
p(f, t) of finding a pair of reactants at the separation f at time t. As shown in Appendix [Bj 
the lattice model considered in this paper leads to 





dt 1 



p(r,t) = V-D B (l-c v (f,t))Vp(r,t)+p(r,t)Vc v (f,t) -k(r)p(r,t) 



(51) 



in the continuous limit, where the diffusion constant is defined by Db = b 2/ -/B and the concen- 
tration of inert particles is denoted by c v (f,t). The first term on the right-hand side includes 



the diffusion term influenced by the concentration gradient of inert particles 
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and the 
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drift term is induced by the spurious potential defined by, 



U = -Jfc B Tln[l - c v (r,t)] 



(52) 



Equation f[5"Tj) can be rewritten in terms of the potential as 



0_ 

Of 



■p(r,t) = V ■D B (l-c v (f,t)) 



Vp (r, t) + p (r, t) 



W 



— k(r)p (r, i) 



(53) 



The perfectly reflecting boundary condition at R is imposed to express that the reactants cannot 
penetrate each other. 

We calculate the escape probability on the basis of Eq. (|53|) . We consider the case that the 
density of the inert particles is stationary and inhomogeneous. Both the intrinsic reaction rate 
and the density of the inert particles are assumed to be isotropic. The equation for the survival 
probability is obtained from Eq. ( |53i) by introducing the adjoint operator as, [39M41] 

(r, t) = V • D B (1 - c v (r))VS (r, t) - D B \vc v (r)l • VS (r, t) - k (r) S (r, t) , (54) 
at L J 

where V in the square brackets operates only on c v (r). The perfectly reflecting boundary 
condition at r = R is represented by, 



0. 



(55) 



42| the 



When the reaction takes place at the reaction radius, R, with the intrinsic rate, ko 
equation for the escape probability defined by tp(r) = lim^oo S (r, t) in <i-dimension satisfies 

U \ d 



^expf-^] ^-D B (l-c v (r))r d 1 exp 



-<p(r) = 0, 



(56) 



k B T J dr \ k B T J dr 

using the potential defined by Eq. ( |52|) . The boundary conditions are given by lim^oo ip(r) = 1 
and 



S d D B {l-c v {R))^(r) 
or 



k (p{R), 



(57) 



r=R 



where the surface area of the d- dimensional sphere is given by Sd = dn d ^ 2 /T((d/2) + 1). For 2 
and 3 dimensions, S2 = 2ttR and S3 = 4nR 2 . The solution of Eq. ( 1561) subject to the above 
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mentioned boundary conditions is obtained as, 

1 Sa 



JR D B (1 - c v Wrt 1 h (1 - c v (R))R d - 1 

¥>(0 = -foo 1 1 • (5J 



R D B {1 - c v (rx)) 2 ^- 1 k (1 - c v (R))R d -i 
In the limit of perfectly absorbing boundary condition, the escape probability simplifies into, 

i , r° , i 



' P(r) = L "'Ul-c^Wrr' L dn (l-Mrt- (59) 

According to Eq. (159)) . the escape probability is independent of the density of the inert 
particles when the density is homogeneous. This result is consistent with that obtained in 
the lattice system in Sec. IV. On the other hand, when the density of the inert particles is 
inhomogeneous, the escape probability in general depends on the density of the inert particles. 
The escape probability is lower than that for the homogeneous density of inert particles if c v (r) 
has a positive slope. The recombination reaction can be assisted by a positive density gradient 
of the inert particles. On the other hand, the recombination can be hindered by a negative 
density gradient of the inert particles. 

VI. CONCLUSIONS 



In this paper, the time evolution equations for the survival probability of a geminate pair 
under the presence of many inert particles are derived and the results are compared to the 
simulation results. If we ignore correlations higher than two-point correlations, Eq. ( Tl9l) is 
derived. In this lowest order approximation, the influence of inert particles is described by using 
the mean field expression of the tracer-diffusion constant in the reaction-diffusion equation. 

In the lowest order approximation, the so-called Bardeen-Herring back correlation is not 
taken into account. The Bardeen-Herring back correlation is the tendency of the preferred 
jump of the diffusing particle back to the previously occupied empty site. We have shown that 
the reaction interferes with the Bardeen-Herring back correlation. In the react ion- diffusion 
equation, the transition operator describing diffusion is influenced by the reaction strength while 
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leaving the reaction term unaltered. If the interference between the reaction and the Bardeen- 
Herring back correlation is taken into account, the reaction-diffusion equation becomes very 
complicated and cannot be solved analytically. By taking into account the Bardeen-Herring 
back correlation but ignoring the interference between reaction and the Bardeen-Herring back 
correlation, we obtain the reaction-diffusion equation given in terms of the improved expression 
of the tracer-diffusion constant. The influence of the excluded volume interactions is taken 
into account solely by the tracer-diffusion constant. The tracer-diffusion constant decreases by 
increasing the concentration of the inert particles since the diffusive motion of reactive species 
is hindered by the presence of the inert particles. 

By comparison of the theoretical results with the results of numerical simulations the in- 
terference between reaction and the Bardeen-Herring back correlation is shown to be small as 
long as overall reaction yield is concerned. The interference also influences the transient decay 
of the survival probability of a geminate pair in the diffusion-controlled limit. When the initial 
concentration of the inert particles is high, we show that interference between reaction and the 
Bardeen-Herring back correlation accelerates the transient decay of the survival probability in 
the diffusion-controlled limit. 

Recently, Schmit et al. studied the reaction in microfluid by the Monte-Carlo simulation. 



261 ] The simulation results are well approximated by the survival probability of a pair of reac- 



tants obtained from the react ion- diffusion equation similar to Eq. (1191) : it is suggested to use 
the reaction-diffusion equation in which the mutual diffusion coefficient is substituted by the 



effective tracer-diffusion coefficient of a tagged particle in a sea of inert particles. [26] They 
calculated the mean first passage time instead of the survival probability. Since the mean first 
passage time can be given in terms of the inverse of the overall reaction rate, their results 
are consistent with our conclusion, although the slightly different expression of the effective 
tracer-diffusion coefficient is used in their equation. According to our work, the approximation 
would fail in the transient decay of the survival probability in the diffusion controlled limit at 
high concentration of inert particles. 

The above conclusions are obtained by assuming the homogeneous distribution of the inert 
particles. We have also formulated a way to obtain the survival probability of a geminate pair 
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when the initial distribution of inert particles is inhomogeneous. The reaction yield is increased 
when the reaction proceeds in the presence of a positive density gradient of the inert particles 
which inhibits the escape of reactants. The effect can be interpreted as a cage effect. Although 
we need further investigation for the kinetics of the survival probability in the presence of the 
density gradient of the inert particles, our results on the escape probability indicate that the 
crowding promotes reactions when the density of the inert particles increases with the distance 
from the location of the immobile reactant. 

Recently, the cage effect of crowding by inert particles was introduced by Kim and Yethiraj 
to interpret the increase of association rate constant with increasing the concentration of inert 
particles obtained by Brownian dynamic simulations when the intrinsic reaction rate constant 
was small. (3] The cage effect in this paper is similar to that introduced by them in the 
sense that the presence of inert particles surrounding a reactant pair promotes the reaction 
but the mechanism is slightly different. In our case, the concentration gradient toward one of 
the reactant pair assists the recombination reaction, while the association rate in their case 
increases by the increase of the recollision probability due to the high concentration of inert 
particles rather than the concentration gradient. 

Appendix A: Diffusion kernel 

By introducing the operator identity, 
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into Eq. (j32p . and from the definition given by Eq. (1521) we derive, 
sG n + b r o \rh», ro, sj — 5 r ^5^ = p)u>G n + b r ° \mm, r + fc a o, s 

+7b(1 — c)G ^f?», n + b r o \ fh — b a »,r + b a o, sj — 74 G n + b r o \m», ro, s 
— pujG n + b r o \mm, b a o : sj + 7 b(1 — c)G (n», n + b r o \ fh — b a », b a o y s 



+ ^2 Sffia \1bcG [ftm, n + b r o \ fh + b a », —b a o, sj + ( 7 t - 7sc) G n + 6 r o \mm, b a o, s 

a 

— k (m) G ( n«, n + 6 r o ro, s ) , (A2) 



where 

It = T«» + Tb(1-c). (A3) 

fn the presence of (m), G (^n»,n + b r o |m«,ro, sj does not satisfy the condition of transla- 
tional invariance of fh against n. 

If we ignore k(rh) in Eq. (1A2I) . the translational invariance is satisfied. The equation for 
G^ ($, b r \ f, s j introduced in Eq. (|38|) is given by, 

sG^ (i, b r \r, s) = ^ [ 7 ^ (T) (7, 6 r |f + 6«, s) + 7fl (l - c)G^ (l- b a , b r \r + b a , s) 

a 

- lt G^ (j,b r \r,s)] -6?^ [fwGW (iX\b a ,s) + lB (l-c)G^ (i-b a ,b r \ba,s] 

a 

+ E % Q [tbcG (t) (f+ 6„, 6 r | - 6 a , s) + ( 7 * - Ibc) (i, b r \b a , s)] • (A4) 



By applying Fourier transformation, 

# (X b r \h, s) = ^ XI exp ^k-i + h- r) G (T) 6 r |r, sj 

G (fc, 6 r |f, s) = X exp |i£ • f G (T) 6 r |r, s' 



(A5) 
(A6) 



we obtain, 

sg ( k, b r | h, s ) = exp ( ih ■ b r ) + \u>t {ba) exp (—ih-b c 



It 



g[k,b r \h,s 



[ ex p (-^ 



z/i • b a ) oj b ( &a ) c + 7s exp ( ih ■ b a ) — u t ( b 



G [k,b r \b a ,s) , (A7) 
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where we define, 



1s=lt- 7bc, 

wb (f) = 7^ exp (ik ■ rj , and u t (f) = 7™ + 7s(l - c) exp ( zTc i ■ f ) . 
It is convenient to introduce the function, 



(A8) 
(A9) 



Q(k,f,s) 



(2tt) c 



tt exp ( — ih ■ r 
d d h 



s + T t -n t (k,h) 



(A10) 



where T t = Idrft and Qt{k, h) is defined by, 



n t (k, h) = 2d lw \{h) + 7B (1 - c)A(fc - fr) 



(All) 



A closed set of equations can be obtained from Eq. (1A7j) as, 
G f fc, fr r sj = b q - b r , s) + ^ \lsQ(k, b q - b a , s) + u B (b^j cQ(k, b q + b a , s)- 

a 

ut (ba) Q(k, b q , s)] G (k, b r \b a , s) . (A12) 



The solution is independent of the position of the reactive sink in this approximation. By 
introducing the solution of Eq. flA12j) into Eq. ( HT|) and using Eq. fl39l) and Eq. (HOI) the 
survival probability is obtained after numerical inverse Laplace transformation and inverse 
Fourier transformation. 

An equation for G^ a ' ( b r |r, s ) is obtained by applying the operator identity, 



s-L 



as 



Hio| 



rwO 



- I 1 H = — L 

s — L rw0 



(A13) 



sG {a) (b r \f, s)-% = [G (a) (X |f + 6 a , s) - GW (& r |f, s)] -^J^G^ (X |&«, s) + 

E [tbcGW (X I - b a , s) + ( 7t - 7sc) G<°> (X |6 a , s)] . (A14) 

The solution depends on a position of an inert particle through its relative vector against the 
initial position of a mobile reactant. 



27 



Appendix B: Excluded volume interactions under inhomogeneous distribution of inert 
particles 



Since the particle jumps to a vacant site, p(f, t) obeys, 

-p(r, t) = 7s (p (f+ bj», f(f), tj —p (r»,r + bj(J), t 



d_ 



(Bl) 



where p \ ?»,r'(f),t\ denotes the joint probability at time t that 
reactive particle and the site r' is empty. We note the relation, [34 



he s ite r is occupied by a 
35, 



43] 



p ( r», r'(j), t ) = p(f, t) — p ( r», r'o, t J — p ( r«, r'», t 



(B2) 

where p (^f»,r'o,tj and p ^r«,r'«,tj denote the joint probabilities at time t that the site r is 
occupied by a reactive particle and the site r' is occupied by an inert particle and that both 
sites are occupied by reactive particles, respectively. Equation (1B1 [) can be rewritten as 

—p (f, t) = 7 B ^2 \P (^*+ bj 9 , tj —p (r«, t) — p (f + bj», fo, t j + p (?•, r + 6jO, t 



— k{r)p (f, t) 



(B3) 



By assuming that the pair correlation function depends on the distance vector between the 



reactant and the inert particle alone, the joint probability function can be written as, 



p ( r», r'o, t ) = p (f,t) c (r',t) a(r — r') 



44 1 



(B4) 

where the occupation probability by an inert particle is denoted by c(r, t). When the spatial 
variation of <r(r) is smaller than p(f, t) and c(f,t), we obtain in the limit of small lattice 
spacing, 

^p (r + bj», fo, tj = [p (f, t) c (r, t) + b 2 c (f, t) V 2 p (f, t)] a(0), (B5) 

j 

^p (f;f+ bjO, tj = [p (f, t) c (r, t) + b 2 p (r, t) V 2 c (f, t)] cx(0), (B6) 
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where lim^o^') = a (®)- By introducing Eqs. ( 1B5I) - (1B6I) . Eq. (1B3I) in the limit of b — > 

becomes, 



a(0)c v (r,t))Vp(r,t) + a(0)p(r,t)Vc v (r,t) -k{r)p{r,t), (B7) 



where D B = b 2 ^ B and the concentration of inert particles is denoted by c v (f, t) in the continuum 
limit. In the mean field approximation in which <r(0) = 1, the derivation follows from that given 



previously. 
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35] Equation (1B7|) can be expressed in terms the correlation factor as, 



Of 



p(r,t)=V-D B {l-c v (r,t))f c Vp(r,t) + a(0)p(r,t)Vc v (r,t) -k(r)p(r,t), (B8) 



where <r(0) = f c + (1 — f c )/c v (r,t). f c is given by Eq. (|50|) . though, strictly speaking, the 
calculation of the correlation factor f c is restricted to the homogeneous distribution of inert 
particles. In the mean field approximation, we have a(0) = 1 and Eq. (1B7I) leads to Eq. (l5"Tj) . 
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